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ABSTRACT 

Time evolution of the expectation values of various dynamical operators of the har- 
monic oscillator with dissipation is analitically obtained within the framework of the Lind- 
blad's theory for open quantum systems. We deduce the density matrix of the damped 
harmonic oscillator from the solution of the Fokker-Planck equation for the coherent state 
representation, obtained from the master equation for the density operator. The Fokker- 
Planck equation for the Wigner distribution function is also solved by using the Wang- 
Uhlenbeck method of transforming it into a linearized partial differential equation for the 
Wigner function, subject to either the Gaussian type or the 5-function type of initial condi- 
tions. The Wigner functions which we obtain are two-dimensional Gaussians with different 
widths. 



1. Introduction 

In the last two decades, the problem of dissipation in quantum mechanics, i.e. the 
consistent description of open quantum systems, was investigated by various authors [1-6]. 
Because dissipative processes imply irreversibility and, therefore, a preferred direction in 
time, it is generally thought that quantum dynamical semigroups are the basic tools to 
introduce dissipation in quantum mechanics. The most general form of the generators of 
such semigroups was given by Lindblad [7-9]. This formalism has been studied for the 
case of damped harmonic oscillators [8,10,11] and applied to various physical phenomena, 
for instance, the damping of collective modes in deep inelastic collisions in nuclear physics 
[12-14] and the interaction of a two-level atom with the electromagnetic field [15]. 

This paper, also dealing with the damping of the harmonic oscillator within the Lind- 
blad theory for open quantum systems, is concerned with the time evolution of various 
dynamical operators involved in the master and Fokker-Planck equations, in particular 
with the time development of the density matrix. In [14] , by using the characteristic func- 
tion of the density operator, we obtained analytical expressions for the first two moments 
of coordinate and momentum of the harmonic oscillator with dissipation. In [16] the Lind- 
blad master equation was transformed into Fokker-Planck equations for quasiprobability 
distributions and a comparative study was made for the Glauber P, antinormal ordering 
Q and Wigner W representations. In [17] the density matrix of the damped harmonic 
oscillator was represented by a generating function, which was then obtained by solving a 
time-dependent linear partial differential equation derived from the master equation. 

The aim of this work is to explore the physical aspects of the Fokker-Planck equation 
which is the c-number equivalent equation to the Lindblad master equation. Generally the 
master equation for the density operator gains considerably in clarity if it is represented 
in terms of the Wigner distribution function which satisfies the Fokker-Planck equation. 
It is worth mentioning that the master and Fokker-Planck equations agree in form with 
the corresponding equations formulated in quantum optics [18-25]. 

The content of this paper is arranged as follows. In Sec. 2 we review briefly the 
derivation of the master equation for the density operator of the harmonic oscillator. In 
order to get an insight into physical meanings of the master equation, we first split it up into 
several equations satisfied by the expectation values of the dynamical operators involved in 
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the master equation. These equations are then solved analytically. In Sec. 3 we transform 
the master equation into the Fokker-Planck equation by means of the well-known methods 
[4,26-30]. We show how the density matrix, subject to the Poisson distribution type 
initial condition, is extracted with the help of the solution of the Fokker-Planck equation 
for the coherent state representation. Then the Fokker-Planck equation for the Wigner 
distribution is solved by the Wang-Uhlenbeck method of transforming it into a linearized 
partial differential equation for the Wigner function, subject to either the Gaussian type 
or the 5-function type of initial conditions. Finally, conclusions are given in Sec. 4. 
2. The master equation for the damped quantum harmonic oscillator 
The rigorous formulation for introducing the dissipation into a quantum mechanical 
system is that of quantum dynamical semigroups [2,3,7]. According to the axiomatic theory 
of Lindblad [7,9], the usual von Neumann-Liouville equation ruling the time evolution of 
closed quantum systems is replaced in the case of open systems by the following equation 
for the density operator p: 

^ = mm (2.D 

Here, $t denotes the dynamical semigroup describing the irreversible time evolution of the 
open system in the Schrodinger representation and L the infinitesimal generator of the 
dynamical semigroup $ t . Using the structural theorem of Lindblad [7] which gives the 
most general form of the bounded, completely dissipative Liouville operator L, we obtain 
the explicit form of the most general time-homogeneous quantum mechanical Markovian 
master equation: 

where 

L(p(t)) = ~[H,p(t)} + ±:52([V jP (t),Vf] + [V^pW;]). (2.3) 
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Here H is the Hamiltonian of the system. The operators Vj and are bounded operators 
on the Hilbert space H of the Hamiltonian. 

We should like to mention that the Markovian master equations found in the literature 
are of this form after some rearrangement of terms, even for unbounded Liouville operators. 
In this connection we assume that the general form of the master equation given by (2.2), 
(2.3) is also valid for unbounded Liouville operators. 



In this paper we impose a simple condition to the operators H, V}, that they 
are functions of the basic observables q and p of the one-dimensional quantum mechanical 
system (with [q, p] = i%) of such kind that the obtained model is exactly solvable. A precise 
version for this last condition is that linear spaces spanned by first degree (respectively 
second degree) noncommutative polynomials in p and q are invariant to the action of the 
completely dissipative mapping L. This condition implies [8] that Vj are at most first 
degree polynomials in p and q and H is at most a second degree polynomial in p and q. 
Then the harmonic oscillator Hamiltonian H is chosen of the form 

P- / n 1 9 mo; 2 9 , n lS 

H = H + ^(pq + qp), H Q = —p 2 + —q 2 . (2.4) 

With these choices the Markovian master equation can be written [10]: 

dp % % % 

jj: = —fi[Ho,p] - 7^(A + fj)[q, PP + pp\ + -^-( A ~^)[P' PQ + QPl 

Dpp -[Q, hP)} ~ [P. \P,P]] + ^([9. \P,P]\ + [?>P]D. (2-5) 



h 2 L "' L "' PJJ h 2 

where D pp ,D qq and D p9 are the diffusion coefficients and A the friction constant. They 
satisfy the following fundamental constraints [10]: 

i) D pp > 0, ii) D qq > 0, iii) D pp D qq - D 2 pq > X 2 h 2 /4. (2.6) 

The equality fx = A is a necessary and sufficient condition for L to be translation invariant 
[8]: [p, L{p)] = L([p, p]). In the following general values for A and \i will be considered. In 
the particular case when the asymptotic state is a Gibbs state 

Pg(oo) = e"^/Tre"^, (2.7) 

these coefficients reduce to 

\ + n hu A — fx h hu 

D rm = Timu) coth - — , D aa = coth - — , D va = 0, (2.8) 

pp 2 2kT' qq 2 mu 2kT' pq K J 

where T is the temperature of the thermal bath. 

Introducing the annihilation and creation operators via the relations 



2mu 



(a+ + a), p = iJ-—(a + -a), (2.9) 



we have H = %u>{a + a + \) and the master equation (2.5) takes the form 

^ = ^(D 1 + p){a + a + p - a + pa+) + ^(D 1 - p)(pa+a + - a + pa+) 

+ \{D 2 + A + iu)(apa + - a + ap) + \{D 2 - A - iu)(a + pa - paa + ) + H.c, (2.10) 
where 

D 1 = UmuD qq -^ + 2iD pq ), D 2 = \{muD qq + ^E). (2.11) 

Ti muj n mu 

In the literature, master equations of the type (2.5) or (2.10) are encountered in 
concrete theoretical models for the description of different physical phenomena in quantum 
optics [18-24], in treatments of the damping of collective modes in deep inelastic collisions 
of heavy ions [12-14,31-34] or in the quantum mechanical description of the dissipation for 
the one-dimensional harmonic oscillator [4,6,8,10,28,29]. A classification of these equations, 
whether they satisfy or not the fundamental constraints (2.6), was given in [17]. 

The meaning of the master equation becomes clear when we transform it into equa- 
tions satisfied by various expectation values of operators involved in the master equation, 

< A >= Tr[p(t)A], where A is an operator composed of the creation and annihilation 
operators. Multiplying both sides of (2.10) by a and taking throughout the trace, we get 

Y < a >= -(A + iu) < a > +p < a + > . (2.12) 

Similarly, the equation for < a + > is given by 

Y < a + >= -(A - iu) < a + > +p < a > . (2.13) 

In the absence of the second term on the right-hand side of (2.12) and (2.13), the two 
equations represent independently a simple equation of oscillation with damping. By 
coupling (2.12) to (2.13) we get a second order differential equation for either < a > or 

< a + >. For example, we obtain: 

d 2 d 

< a > +2A— < a > +(A 2 + uj 2 - p 2 ) < a >= 0, (2.14) 

which is the equation of motion for Brownian motion of a classical oscillator, but without 
the term corresponding to random process. Because of the vanishing terms on the right- 
hand side, we may equally state that (2.14) is the equation of motion with zero expectation 



value of the random process. In the study of the Brownian motion of a classical oscillator, 
one replaces the second-order differential equation of motion of the type (2.14) by two 
equivalent first-order differential equations [35] which are precisely the Langevin equations. 
Accordingly, we may say that (2.12) and (2.13) are the Langevin equations corresponding to 
(2.14), but without the random process term. The integration of (2.14) is straightforward. 
There are two cases: a) fx > u> (overdamped) and b)/i < to (under damped). In the case a) 
with the notation v 2 = [i 2 — cu 2 , we obtain: 

< a (t) >= e~ xt [< a(0) > (coshz/t - i- sinhz4) + - < a+(0) > sinhi/t]. (2.15a) 

In the case b) with the notation Q? = u 2 — fx 2 , we obtain: 

< a(t) >= e~ xt [< a(0) > (cosfit - sin fit) + ^ < a + (0) > sin fit]. (2.15b) 

Here < a(0) > or < a + (0) > stands for the initial value. The expression for < a + (t) > 
can be obtained simply by taking the complex conjugate of the right-hand side of (2.15). 

For the computation of quantal fluctuations of the coordinate and momentum of the 
harmonic oscillator, we need the expectation values of product operators, such as a +2 , a 2 or 
a + a. The dynamical behaviour of these product operators can be well surveyed by deriving 
the equations satisfied by their expectation values. By following the same procedure as 
before, employed in the derivation of (2.12), we find: 

d 

— < a 2 > +2(A + iu) < a 2 >= 2/j < a + a > +D 1 + fi, (2.16) 

(Jib 

< a+ 2 > +2(A - iu) < a +2 >= 2^ < a+a > +D{ + fi, (2.17) 

< a + a > +2A < a + a >= //(< a +2 > + < a 2 >) + D 2 - A. (2.18) 

Eqs.(2.16)-(2.18) exibit clearly that the expectation value of one double operator has effects 
upon the dynamical behaviour of other double operators due to the presence of non-zero 
terms on the right-hand sides. The solutions of these equations are readily obtained by 
transforming them into two differential equations satisfied by the sum and the difference 
of two double operators, < a 2 > and < a +2 >. We obtain: 

^(< a 2 > + < a+ 2 >)+4\j t (< a 2 > + < a+ 2 >) +4(A 2 + u 2 -fi 2 )(< a 2 > + < a+ 2 >) 



= \ [(A + fi)muD qq -(\- fl )£™+ 2ivD pq ] = D, (2.19) 

±(<a 2 >-< a+ 2 >) + 2A(< a 2 > - < a+ 2 >) + 2iuj(< a 2 > + < a+ 2 >) = 4i^. 

(2.20) 

The solution of (2.19) is straightforward and with the help of which both equations (2.18) 
and (2.20) can be immediately solved. We find: 

< a 2 >= e"^[(l -i-)C ie ^ + (l+i^)C 2 e- 2 ^ - z^C 3 ] + £gf % + »5f , (2.21a) 

v v uj 8A(A / — v z ) nX 

< a+a >= e-^fyc^" - C 2 e~ 2 ^) + C 3 ] + + (2.22a) 

for the overdamped case fx > uj and 

< o2 > =e - [ (C 1+i C^)cos2m + (C 2 - i C^)sin2m-i5c3] + i ^| y+i ^, 

(2.21b) 

< a+a >= e- 2At [^(C! sin20t - C 2 cos20t) + C 3 ] + 8A(A f^ Q2) + (2 ' 22b) 

for the underdamped case uj > fi. The expression for < a +2 > can be obtained by 
taking the complex conjugates of (2.21a) and (2.21b). Here, Ci,C 2 ,Cs are the integral 
constants depending on the initial expectation values of the operators under consideration. 
In particular, if D qq = D pq = and fi = A we obtain for the underdamped case uj > fx 
the same equations as those written by Jang in [31] for his model on nuclear dynamics 
based on the second RPA at finite temperature. If the condition fx <C uj is satisfied 
throughout in the underdamped case, the expectation values (2.21b) and (2.22b) reduce 
to (Ci + iC 2 ) exp[-2(A + iuj)t] + 2(A f_^ 2) (A - iu) and C 3 exp(-2At) + ^§^, respectively. 
These are the underdamped solutions of (2.16) and (2.18), but in the right-hand sides 
setting fx = 0. For time t — > 00, we see from (2.22) that 

<° + °>= 8A(*H%-,.) + H^ < 2 ' 23 > 

In the particular case when the asymptotic state is a Gibbs state (2.7), (2.8) and fi = A, 
we get 

< a+a >= -(coth — - 1) = (exp — - l)" 1 =< n >, (2.24) 
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which is the Bose distribution. This means that the expectation value of the number 
operator goes to the average thermal-phonon number at infinity of time. From the identity 

oo 

< a + a >= ^ m < m\p(t)\m > (2.25) 

m=0 

it follows 

< n > m 

< m\p(oo)\m>= r—r. r . (2.26) 

'' v (1+ < n >) m+1 v ' 

In deriving this formula, we have made use of the identity J2m=o mxm = x/(l — x) 2 . The 

expression (2.26) shows that in the considered particular case the density matrix reaches its 

thermal equilibrium - the Bose-Einstein distribution, whatever the initial distribution of 

the density matrix may be. When the initial density matrix < m\p(0)\m > is represented 

by a distribution of the form N m /(1 + iV) m+1 , where N stands for the average phonon 

number, the relation (2.25) implies that < a + a >= N. When the initial density matrix is 

characterized by a distribution of the form 

ml 

(2.25) implies that < a + a > becomes also N. Eq.(2.27) is nothing but a Poisson distri- 
bution. If the initial density matrix is represented by a Kronecker delta 5 ms , we see from 
(2.25) that < a + a >= s, which corresponds to the initial s-phonon state. 

For the master equation (2.10) of the harmonic oscillator, physical observables can 
be obtained from the expectation values of polynomials of the annihilation and creation 
operators. So, for the position and momentum operators q and p via the relations (2.9), 
we can evaluate either the second moments or variances (fluctuations), by making use of 
the results (2. 15), (2. 21), (2. 22). Such quantities have been also calculated [10,14], by using 
other methods. 

3. Fokker-Planck Equations 

One useful way to study the consequences of the master equation (2.10) for the den- 
sity operator of the one-dimensional damped harmonic oscillator is to transform it into 
more familiar forms, such as the equation of motion for the density matrix [17] or the 
equations for the c-number quasiprobability distributions Glauber P, antinormal ordering 
Q and Wigner W associated with the density operator [16]. In the second case the result- 
ing differential equations of the Fokker-Planck type for the distribution functions can be 



solved by standard methods [26-29] employed in quantum optics and observables directly 

calculated as correlations of these distribution functions. 

3. 1 Calculation of the density matrix from the Fokker-Planck equation 

The Fokker-Planck equation, obtained from the master equation and satisfied by the 

Wigner distribution function W(a, a*,t), where a is a complex variable, has the form [16]: 

dW(a^a ,t) = _{d_[_^ x + ^ + ^ + A- [-(A - iu)a* + /ia]}W(a, a*, t) 

When we substitute the P representation function P(a,a*,t) for W(a,a*,t) and the co- 
efficients D\ + /j, for D\, D2 — A for D 2 in the above equation, we get the Fokker-Planck 
equation for the coherent representation [16]. Despite the formal similarity to averaging 
with a classical probability distribution, the function P(a,a*,t) is not a true probability 
distribution. Because of the overcompleteness of the coherent states, the P function is not 
a unique, well-behaved positive function for all density operators. Cahill [36] studied the P 
representation for density operators which represent pure states and found a narrow class 
of states for which the P representation exists. They can be generated from a particular 
coherent state \a > by the application of a finite number of creation operators. Also Cahill 
[37] introduced a representation of the density operator of the electromagnetic field that is 
suitable for all density operators and that reduces to the coherent state P representation 
when the latter exists. The representation has no singularities. 

The Fokker-Planck equation (3.1) or the similar equation for the P representation, 
subject to the initial condition 

P(a,a*,0) =5(a-a )5(a* -c$), (3.2) 

where «o is the initial value of a can be solved [30,38,39] and the solution to (3.1) (Green 
function) for the P representation is found to be 

2 1 

P(a,a*,t) = =^=exp{--— -— [cr 22 (a;-ao) 2 +o-ii(a*-ao) 2 - 2cr i2|«-«o| 2 }, (3.3) 

ny / deta(t) aeta-(tj 

where a = (cr^), 

a ij{ t ) = ^2 ftisSjr ~ b is (t)b jr (t)]o- sr (oo). (3.4) 
s,r=l,2 



The function cko and its complex conjugate, which are still functions of time, are given by 



"o = hi(t)a + b 12 (t)aQ. (3.5) 
The functions bij obey the equations 

i=i,2 

with the initial conditions bj s (0) = Sj s and cr(oo) is determined by 

Ca(oo) + a(oo)C T = Q p , (3.7) 

where 

C=( X + lU V Q P =fn 1+ x n*7 A V ( 3 - 8 ) 

We get 

6 n =6^ 2 = e _A *(cosOt-z^sinOt), b 12 = b 21 = ^e~ xt sinQt, (3.9) 

with O 2 = u 2 — p 2 . While the functions <Tn,cr 22 and ao are complex with an = a^ 2 , the 
functions deter (t) and <J\ 2 are real. 

The solution of the Fokker-Planck equation has been written down providing the 
diffusion matrix Q p is positive definite. However, the diffusion matrix in the Glauber P 
representation is not, in general, positive definite. If the P distribution does not exist as a 
well-behaved function, the so-called generalized P distributions can be taken that are well- 
behaved, normal ordering functions [40] . The generalized P distributions are nondiagonal 
expansions of the density operator in terms of coherent states projection operators. 

In the coherent representation [41,42] the density operator p(t) is expressed by 

p(t) = J P(a,a*,t)\a>< a\d 2 a, (3.10) 

where d 2 a = d(Rea)d(lma) and \a > is the coherent state. The matrix element of p{t) 
in the n quantum number representation is obtained by multiplying (3.10) on the left by 
< to | and on the right by \n > . By making use of the well-known relation 

/ 1 1 i2n or . . . 

\a >= exp(--|a| ) > —j=\ n >> ( 3A1 ) 
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we get 

< m\p(t)\n >= J: I a n a* m P{a,a\t)evp{-\a\ 2 )d 2 a. (3.12) 

ym\n\ J 

Upon introducing the explicit form (3.3) for P(a,a*,t) into (3.12), we obtain the desired 
density matrix for the initial coherent state. However, due to the powers of complex 
variables a and a* in the integrand, the practical evaluation of the integral in (3.12) is not 
an easy task. Instead, we use the method of generating function [17,31] which allows us to 
transform (3.12) into a multiple-differential form. When we define a generating function 
F(x,y,t) by the integral 

F(x,y,t) = J P{ot, a*, t) exp(— \a\ 2 + xa + ya*)d 2 a, (3.13) 

we see that the density matrix is related to the generating function by 

< m\p(t)\n >= -^={^r{^) n F{x^t)\ x=y=Q . (3.14) 
Vmini ox ay 

Since the P representation is in a Gaussian form, the right-hand side of (3.13) can be 
evaluated analitically by making use of the identity 

/ exp(-aH 2 + bz + cz* + ez 2 + fz* 2 )d 2 z = , = exp ahc + h2 f +f e (315) 
J Va 2 - Aef ' a 2 -4ef 

which is convergent for Re a > |e* + /|, while 6, c may be arbitrary. We find: 
2 1 

F = —= exp{xy - — [a u (x - o^) 2 + a 22 (y - «o) 2 ~2(a 12 + 2){x - a* Q )(y - a )]}, (3.16) 
\/L4 A 



where A = d — A{a\ 2 + 1), d = deter = o\\o 22 — o\ 2 . A formula for the density matrix can 
be written down by applying the relation (3.14) to the generating function (3.16). We get 



< m\p{t)\n >= 2 __ exp{[- -[(T22Q0 + ff ii"o 2 ~ 2 (^i2 + 2)|a | 2 ]} 
■s/\A\ A 

( _ 1 ^ n\ +n 2 nm+n— 2(m +n 2 +n 3 ) 

x Y i r ,7 " ttt (3.17) 

ni\n 2 \ns\(m — 2ri\ — n^yAn — 2n 2 — n^y. 

"1,12,13=0 

where 

a^a%(d-2a 12 ) n *[a 11 a* - (a 12 + 2)« ] m - 2ni " n3 [a 22 a - 2(<r 12 + 2)o^]— 2 "™ 



E = 



j^m+n-(n 1 +n 2 +n 3 ) 
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The expression (3.17) is the density matrix corresponding to the initial coherent state. At 
time t = 0, the functions an, 022 and a 12 vanish and «o goes to a®. In this case the density 
matrix reduces to 

< m\p(0)\n >= «gftQ m exp(-|a | 2 ), (3.18) 

\Jm\n\ 

which is the initial Glauber packet. For the diagonal case the initial density matrix becomes 
the Poisson distribution. At infinity of time, the density matrix (3.17) goes to the Bose- 
Einstein distribution 

< n > m 

< m\p(oo)\n >= - - >)m+1 <W (3.19) 

For other specific initial conditions, it is more convenient to proceed to solve directly the 
master equation in order to extract a closed form for the density matrix [17]. 

3.2 Wigner distribution function 

The Fokker-Planck equation (3.1) can also be written in terms of real coordinates x\ 
and X2 (or the averaged position and momentum coordinates of the harmonic oscillator) 
defined by a = x\ + 1x2, a* = x\ — 1x2, as follows [16]: 

dW . d 1 „w d 



= E A v£r( x i w )+ 1 - J2 QYij^r-W, (3-20) 



dt ^ J dxi y J ' 2 dxidxi 



where 

Since the drift coefficients are linear in the variables x\ and X2 and the diffusion coefficients 
are constant with respect to x\ and X2, (3.20) describes an Ornstein-Uhlenbeck process 
[35,43]. Following the method developed by Wang and Uhlenbeck [35], we shall solve this 
Fokker-Planck equation, subject to either the wave-packet type or the 5-function type of 
initial conditions. By changing the variables x\ and X2 of (3.20) via the relations 

z\ = ax\ + 6x2, z 2 = cx\ + dx2, (3.22) 

the Fokker-Planck equation (3.20) is transformed into the standard linearized partial dif- 
ferential equation [35] expressed as 



dW(zi,Z2,t) . d d 1 ^ d 2 > . 

lH = <" W - + 2 .. Z^g^m^ *.«)■ (3-23) 



M=l,2 
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In deriving this equation we have put a = c* = M J Q , b = d = = v 2 = —A — zfi and 
£>n = D 2 * 2 = ^[(/x - iQ) 2 mD qq + 2(/x - + 

D 12 = £>2i = ^M w + — ^ + (3-24) 
ft ^ u} muj 

1) When the Fokker-Planck equation for the coherent state representation is subject 

to the initial condition 8(a — cto)8(a* — ckq), then the use of the relation between the Wigner 

distribution function and P representation 

W(a, a*,t) = ^ J P(p, p*, t) exp(-2|« - (3\ 2 )d 2 (3 (3.25) 

leads to a Gaussian form for the initial Wigner distribution function at t = 0. The initial 
Wigner distribution function can in turn be related to the Wigner distribution function 
expressed in terms of xio and x 2 o - the initial values of x\ and x 2 at t = 0, respectively. 
Then we have the relation 

2 

W w (x u x 2 , 0) = W w (a, a\ 0) = - exp(-2|a - a \ 2 ). (3.26) 

The right-hand side of (3.26) can be expressed in terms of real coordinates x\ and X2 and 
we thus get the expression which corresponds to the initial condition of a wave packet [41]. 
Accordingly, we now look for the solution of the Fokker-Planck equation (3.20) subject to 
the initial condition 

W w ( Xl , x 2 , 0) = - exp{-2[(x! - x 10 ) 2 + (x 2 - x 20 ) 2 ]}. (3.27) 

By changing the variables x\ and x 2 into z\ and z 2 , {z\ = z\ = z), this initial condition is 
seen to be transformed into 

2 2oj 2 

W w (z, z\ 0) = - exp{— f [q(z - z ) 2 + q*(z* - z* ) 2 - \z - z \ 2 ]}, (3.28) 

where zq is the initial value of z and q = niEt^R^ The solution of (3.20) subject to the 
initial condition (3.28) is found to be 
90 1 

W w (z, z*,t) = == exp{--— [g 2 (z - z e^) 2 + 9l (z* - z* e^) 2 - g 3 \z - z e^\ 2 }}, 

HUJ^f\B w \ ^-Dw 

(3.29) 
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where 

&w - 9i92- ^93^91 - 92 - 1 e +_ 2^"' e i J'fi , 3-e +— — (1-e j. (3.3UJ 
In terms of real variables xi and £2 we have: 

W w (x!,X2,t) = ^==exp{--^-[(j) w (xi-X 1 ) 2 +ij w (x2-X2) 2 + Xw(xi-X 1 )(x2-X2)]}, 

ny/\B w \ 2B W 

(3.31) 

where 

</>w = 9ia* 2 + 920 2 - (fe, ip w = 9i+92- #3, X™ = 2(#ia* + # 2 a) - #3 (a + a*). (3.32) 
The functions x\ and #2, which are also oscillating functions, are given by 

x\ = e~ xt [xio(cos fit + ^ sin fit) + ^20^ sin fit], 

X2 = e~ xt [x2o(cos fit — ^ sin fit) — ^10^ sin fit]. (3.33) 

2) If the Fokker-Planck equation (3.20) is subject to the 5-function type of initial 
condition, the Wigner distribution function is given by 

2fi 1 

W(z, z\t) = _ exp{- — [f 2 (z - z e^) 2 + h{z* - z* e^) 2 - 2f 3 \z - z e^\ 2 }}, 

nu)^/\B\ 

(3.34) 

where 

B = hh~fl /i = / 2 * = ^(e 2 ^-l), /3 = |f(l-e- 2 ^). (3.35) 
In terms of real variables x\ and X2 we have: 

W(x 1 ,x 2 ,t) = \= exp{-^—[(p d (xi - xi) 2 + i) d (x2 - x 2 ) 2 + Xd(xi - xi)(x 2 - x 2 )]}, 

2iry/\B\ Zjd 

(3.36) 

where 

0d = M* 2 + f 2 a 2 - f 3 , if> d = fi + f 2 - h, Xd = 2(/ia* + f 2 a) -f 3 (a + a*). (3.37) 

So, one gets a 2-dimensional Gaussian distribution with the average values x\ and X2 and 
the variances 4>d^d and Xd- Obviously, at time t = 0, all functions £?, fx, $2, fs vanish and 
the Wigner distribution function goes to the 5-function. 
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When time t goes to infinity, x\ and vanish and we obtain the steady state solution: 

W(x u x 2 )= 1 exp[-l ^ W )^(oc)x iXj ]. (3.38) 

27r v /deta M/ (cx)) 2 

The stationary covariance matrix ^(oo) can be determined from the algebraic equation 

Aa w (oo) + a w (oo)A T = Q w . (3.39) 



We obtain: 

W( ^ _ ( 2A (A + V) + ^ 2 )Qn + ^22 + 2u;(A + ^)QS^ 
an " 4A(A 2 + u; 2 -^ 2 ) ' 

(00) = 4A(A 2 + ^-^) ' (3 - 40) 

w -u(X + /Qgg + c^(A - n)Qf 2 + 2(A 2 - ii 2 )QY 2 
ai2(00) = 4A(A 2 + c 2 -^) • 

4. Conclusions 

The Lindblad theory provides a selfconsistent treatment of damping as a possible 
extension of quantum mechanics to open systems. In the present paper we have studied the 
one-dimensional harmonic oscillator with dissipation within the framework of this theory. 
We have carried out an extended calculation of the expectation values of various dynamical 
operators involved in the master equation, especially the first two moments and the density 
matrix. Generally, the time evolution of the density matrix as well as the expectation 
values of dynamical operators are characterized by complex functions with an oscillating 
element exp(±iA/co> 2 — fi 2 t) multiplied by the damping factor exp(— At). We deduced the 
density matrix from the solution of the Fokker-Planck equation for the coherent state 
representation, obtained from the master equation for the density operator. For a thermal 
bath, when the asymptotic state is a Gibbs state, a Bose-Einstein distribution results as 
density matrix. Also for the case that the initial density matrix is chosen as a Glauber 
packet, a simple analytical expression for the density matrix has been derived. The density 
matrix can be used in various physical applications where a Bosonic degree of freedom 
moving in a harmonic oscillator potential is damped. For example, one needs to determine 
nondiagonal transition elements of the density matrix, if an oscillator is perturbed by a 
weak electromagnetic field in addition to its coupling to a heat bath. From the master 
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equation of the damped quantum oscillator we have also derived the corresponding Fokker- 
Planck equation in the Wigner W representation. The Fokker-Planck equation which 
we obtained describes an Ornstein-Uhlenbeck process. We have solved this equation by 
using the Wang-Uhlenbeck method of transforming it into a linearized partial differential 
equation for the Wigner function, subject to either the Gaussian type or the 5-function type 
of initial conditions and showed that the Wigner functions are two-dimensional Gaussians 
with different widths. 
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